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ABSTRACT 

We present time dependent modeling based on the accretion disk limit cycle 
model for a 270 d light curve of the short period SU UMa-type dwarf nova V344 
Lyr taken by Kepler. The unprecedented precision and cadence (1 minute) far 
surpass that generally available for long term light curves. The data encompass 
two superoutbursts and 17 normal (i.e., short) outbursts. The main decay of 
the superoutbursts is nearly perfectly exponential, decaying at a rate ~12 d 
mag~^, while the much more rapid decays of the normal outbursts exhibit a faster- 
than-exponential shape. Our modeling using the basic accretion disk limit cycle 
can produce the main features of the V344 Lyr light curve, including the peak 
outburst brightness. Nevertheless there are obvious deficiencies in our model light 
curves: (1) The rise times we calculate, both for the normal and superoutbursts, 
are too fast. (2) The superoutbursts are too short. (3) The shoulders on the 
rise to superoutburst have more structure than the shoulder in the observed 
superoutburst and are too slow, comprising about a third to half of the total 
viscous plateau, rather than the ~10% observed. However, one of the cccoid cthot 
interpolation schemes we investigate (one that is physically motivated) does yield 
longer superoutbursts with suitably short, less structured shoulders. 
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1. Introduction 

Dwarf novae constitute a subclass of the cataclysmic variables (Warner 1995b), semi- 
detached interacting binaries in which a Roche-lobe filling secondary transfers matter to a 
more massive white dwarf (WD). The transferred material possesses angular momentum, 
and therefore can only "fall" toward the primary down to a radius determined by its specific 
angular momentum (Lubow & Shu 1975). This angular momentum barrier is overcome by 
angular momentum transport within the accretion disk which carries angular momentum 
outward and mass inward. Current thinking about the physical mechanism responsible for 
accretion centers on the magneto-rotational instability (=MRI; Balbus & Hawley 1998) in 
which the shearing amplification of a weak seed magnetic field leads to turbulent transport. 

V344 Lyr is a short period dwarf nova, of subtype SU UMa (Porb = 2.1 h; Still et al. 
2010) exhibiting a supercycle pattern (109.6 d; Kato et al. 2002) consisting of superoutbursts 
separated by normal outbursts. Superoutbursts are long outbursts with slow, exponential 
decays exhibiting superhumps — modulations in the light curve at a period greater than the 
orbital period by a few percent. The data presented in Kato et al. reveal only the brightest 
portions of the superoutbursts; the fainter portions and all of the normal outbursts lie below 
the lower limit of detection. 

The Kepler data set for V344 Lyrae {Kepler ID 7659570) — 270 d at 1 min cadence 
— provides an extraordinary resource for accretion disk modelers. The accretion disk limit 
cycle mechanism, which is currently employed to account for dwarf nova outbursts, is based 
on the storage of material during quiescence in a non-steady-state configuration, followed 
by a dumping of matter onto the accreting WD when a critical surface density is attained 
in the disk, producing an outburst. During outburst the disk is roughly in steady state, 
except for the action of a cooling front that begins at the outer disk and moves inward. 
The thermal-viscous limit cycle model is not a complete model for the accretion disk, but 
rather primarily a model for the outbursts; for instance it fails to account for the high X-ray 
luminosities seen in dwarf novae in quiescence (Wheatley, Mauche, & Mattel 2003, Mukai, 
Zietsman, & Still 2009). 

Previous attempts at constraining physical parameters of the model have relied on 
AAVSO (American Association of Variable Star Observers) and RASNZ (Royal Astronomical 
Society of New Zealand) data for the long term light curves. These data are typically given 
as Id means, with large attendant error bars (~0.3— 0.5 mag). The Kepler data provide an 
entirely new level of precision, and allow detailed model constraints. The viscosity parameter 
a (Shakura & Sunyaev 1973) sets the disk time scales. To first order, ctcoid controls the 
recurrence time for outbursts, and cchot sets the duration of the outbursts. Previous efforts, 
for instance by Smak (1984=S84), utilized the "Bailey relation" relating the rate of the fast 
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decay in DNe with orbital period to infer that a^ot — 0.1 — 0.2. This remains the only 
firm constraint we have in astronomy for a in ionized gas, and has served as a benchmark 
value for modelers calculating the efficiency of the MRl in accretion disks (e.g., Hirose et al. 
2009). By comparing theoretical modeling of the MRl with the observational constraints. 
King, Pringle, & Livio (2007) showed that theoretical values tend to be at least an order 
of magnitude smaller. Most of the simulations have studied the growth and subsequent 
nonlinear saturation of the MRl using idealized "shearing boxes" without realistic vertical 
structure, and for such calculations the asymptotic strength of the magnetic field decreases 
with increasing resolution (Promang & Papaloizou 2007). Vishniac (2009) argues that current 
models are not sophisticated enough in terms of their dynamic range and treatment of the 
vertical disk structure to allow a realistic treatment of such important effects as magnetic 
buoyancy and vertical transport of magnetic helicity by turbulent eddies, and therefore 
comparisons of theoretical a values with those inferred from observations are specious at the 
present time. 

Our goal in this work is to find model parameters that can account for the gross prop- 
erties of the long term V344 Lyr light curve. There will be additional features in the light 
curve that can be compared to the detailed model light curves in order to gain insight into 
shortcomings of the model. We also investigate systematic effects such as the number of grid 
points and time step size. 

In §2 wc present an overview of the physics of accretion disks. In §3 we discuss the 
V344 Lyr light curve, in §4 wc present a review of viscous decays in several well-studied 
dwarf novae, in §5 we detail the numerical model, in §6 we present the results, in §7 we 
discuss the results in a larger context, and in §8 we sum up. 

2. Accretion Disk Physics 

2.1. Axisymmetric Models 

In the accretion disk limit cycle model, gas accumulates in quiescence and accretes onto 
the central object in outburst (e.g., Cannizzo 1993a, Lasota 2001 for reviews). The phases 
of quiescence and outburst are mediated by the action of heating and cooling fronts that 
traverse the disk and bring about phase transitions between low and high states, consisting 
of neutral and ionized gas, respectively During quiescence, when the surface density S(r) 
at some radius within the disk exceeds a critical value Sinax(^")i a transition to the high 
state is initiated; during outburst, when S(r) drops below a different critical value Smin('"), 
a transition to the low state is initiated. Low-to-high transitions can begin at any radius. 
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whereas high-to-low transitions begin at the outer disk edge. This situation comes about 
because in the outburst disk E(r) oc r~^/^ (roughly), and the critical surface densities both 
increase with radius. Since the disk mass accumulated in quiescence is bounded by Sinax(?') 
and Sniin(r), one can define a maximum disk mass 



which will bound the general, time dependent disk mass. 

For normal, "short" outbursts, only a few percent of the stored gas accretes onto the 
central object: the thermal time scale of a thin disk is short compared to the viscous time 
scale, and the coohng front that is launched from the outer edge of the disk almost as soon as 
the disk enters into outburst traverses the disk and reverts it back to quiescence. For disks 
that have been "filled" to a higher level with respect to Mdisk,max, the surface density in the 
outer disk can significantly exceed the critical surface density Emin- In order for the cooling 
front to begin, however, the outer surface density E (router) must drop below Smin (router)- 
Disks in this state generate much longer outbursts, with slower "viscous" plateaus, because 
the entire disk must remain in its high, completely ionized state until enough mass has 
been lost onto the WD for the condition S (router) < Snim(?"outer) to be satisfied. The slow 
decay associated with a superoutburst is a direct reflection of the viscous time scale in the 
outer disk, whereas the subsequent faster decay reveals the thermal time scale. Although 
the decays of dwarf novae at a given orbital period are fairly uniform, and form the basis 
for the Bailey relation, the rise times show a greater variety, reflecting the fact that the 
outburst can be triggered anywhere in the disk (Cannizzo 1998b). Inside-out bursts tend 
to produce slow-rise times, whereas outside-in bursts prodiicc fast rises (Cannizzo, Wheeler, 
& Polidan 1986). Van Paradijs (1983) undertook a statistical study of outburst properties 
of dwarf novae spanning a range in orbital period, from < 2 hr to ~10 hr. By studying 
separately the behavior of long and short outbursts versus orbital period, he found that 
the apparent strong contrast between superoutbursts and normal outbursts in the SU UMa 
systems is due to the fact that the widths of the short (narrow) outbursts scales with orbital 
period, whereas the widths of the long (wide) outbursts is relatively constant. Van Paradijs' 
argument is strengthened by the more extensive data set examined by Ak, Ozkan, & Mattel 
(2002) in which the linear trend in the narrow width- Porb relation is more obvious (compare 
Fig. 4 from Ak et al. with Fig. 3b from van Paradijs). In these studies the widths of the 
outbursts were "corrected" by subtracting the rise and decay times. 




(1) 



and a minimum disk mass 




(2) 



For any dynamical system in which a variety of timescales contribute to a physical 
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observable, in this case the flux from the accretion disk, the controlling time scale will be 
the slowest one. Since all the disk timescales — dynamical, thermal, and viscous — increase 
with radius, the time scale associated with the viscous plateau provides direct information 
about the viscous time scale tvisc at the outer edge of the disk. 

Warner (1995ab) presents a simple physical model for the state of the disk during 
superoutburst. During the superoutburst the disk is roughly in steady state with S(r) oc 
r~^/^ out to some which represents the radius at which the outward moving heating front 
stalls. The rate of accretion within the disk is much greater than the mass transfer rate from 
the secondary star, which can be neglected during this time. Therefore the mass of the hot, 
accreting disk is 

Mdisk ~ yr^S^i„(rd). (3) 

From the standard accretion disk theory relating effective temperature with accretion rate 
in a steady-state disk (Shakura & Sunyaev 1973), one has 

Mdis^ = -^^T^\ird)rl (4) 

where Tesird) = TesC^mmlrd]) is the effective temperature associated with Smin('"d)- Com- 
bining eqns. [3] and |4] (same as eqns. [15] and [16] from Warner 1995a, and [3.46] and [3.47] 
from Warner 1995b) yields 

Mdisk 3G'MiS^i„(r,) ' ^ ^ 

Since the right hand side (excluding dt) is constant, the solution for the mass of the hot 
disk, and also its derivative, is a decreasing exponential function. Using power law scalings 
for Tes{rii) and Smin('"(i), setting equal to the 3:1 resonance radius with the binary orbital 
period, and noting that luminosity varies as Mdisk, one can show (cf. eqn. [20] of Warner 
1995a) 

t^isc ^ 17 d a~T Pi" m}/^ (6) 

where a_i = a/0.1, Ph is the orbital period in hr, and mi = Mi/Mq. (The a in this context 
is the hot state value.) Within the viscous plateau portion of a superoutburst, the instigation 
and propagation of a cooling front at the outer edge is thwarted by virtue of excessive mass in 
the disk. This phenomenon is not new and has been seen in numerous previous studies (e.g., 
Cannizzo 1993b=C93b, Hameury et al. 1998=H98, Ludwig & Meyer 1998, Buat-Menard, 
Hameury, & Lasota 2001, Schreiber, Hameury, & Lasota 2003). The viscous time does not 
represent the duration of the outburst itself, but rather the time scale (approximately in d 
mag"^) associated with the viscous plateau; it represents an e— folding time for the decrease 
of disk mass. 
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2.2. Non-axisymmetric Models 

Whitehurst (1988) discovered a hydrodynamical instability that can occur in the outer 
parts of accretion disks when the mass ratio q = M2/M1 is less than 1/4. For these extreme 
mass ratio systems, expected in CVs below the period gap if the secondary fills its Roche 
lobe, the outer part of the disk can lie beyond the point of 3:1 resonance with the binary 
period. The outer disk can then be excited by the inner Lindblad resonance (Lubow 199 lab), 
causing the global disk oscillation mode that is observed as common superhumps. Before the 
onset of the superhump oscillation, there are well-known spiral dissipation waves in the disk 
that are fixed in the co-rotating frame. Once the superhump oscillation begins, these spiral 
arms advance 180° every superhump cycl^ As a given arm passes between the two stars, 
it expands outward in the shallower potential, but as fluid in that arm compresses back into 
the disk as the disk is maximally distorted, viscous dissipation causes a peak in the photo- 
metric light curve (Simpson & Wood 1998). A secondary source of the superhump signal 
is caused by the variable distance the mass stream must fall from the LI point to the edge 
of the non-axisymmetric, flexing disk. Because the outer disk is thought to expand during 
superoutburst, potentially beyond the 3:1 radius, the "tidal instability" model provides a 
natural explanation for superhumps. If the mass ratio is extreme enough, superhumps can 
also occur outside of superoutbursts because the disk can remain extended beyond the 3:1 
radius, and therefore continue oscillating even into quiescence. Osaki (1989ab) combined 
the accretion disk limit cycle model with Whitehurst's tidal instability model and proposed 
the "thermal-tidal" instability model (TTIM) for superoutbursts. This model posits that 
the presence of an oscillating precessing disk, triggered by the migration of gas beyond the 
3:1 resonance radius, also leads to a greatly enhanced tidal torque acting on the disk. This 
ultimately contracts the outer disk and leads to a superoutburst. Thus in the TTIM one has 
the usual hysteresis relation between surface density and temperature in the disk, and also 
a hysteresis relation between the total angular momentum of the disk and the tidal torque 
acting on the outer disk. (The latter hysteresis is somewhat speculative.) This model was 
examined in detail and compared to the enhanced mass transfer model (EMTM) for super- 
outbursts by Ichikawa, Hirose, &; Osaki (1993) and later workers, most recently Schreiber et 
al. (2004). 

For completeness, we note that Smak has recently challenged the conventional wisdom 
that superhumps are due to a precessing, eccentric disk (Smak 2009abcd). Smak (2009a) 
presents a re-analysis of data taken by earlier workers on Z Cha, WZ Sge, OY Car, and 
lY UMa and claims that a variety of errors, including miscalculated beat phases and an 



^See, for example, www. astro. fit.edu/wood/visualizations. html. 
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assumption that observed eclipses are pure disk eclipses, led to erroneously large values of 
the eccentricity e. Smak's re-analysis of the disk eccentricity in Z Cha leads to e = 0.05 ± 
0.05. Smak (2009b) compares in detail observed values for the amplitudes of superhumps 
with theoretical values determined from bolometric light curves produced by smooth-particle 
hydrodynamics calculations. He finds the observed values to be ~10 times larger than the 
theoretical values. Smak (2009c) presents a re-analysis of five eclipses of Z Cha observed by 
Warner & O'Donoghue (1988). Smak examines two local minima at orbital phases $ ~ —0.05 
and 0.04. He argues that the first minimum is not due to an occultation but rather absorption 
in the mass overflow stream, and the second minimum coincides with the trajectory of the 
mass overflow stream. The upshot of his analysis is that superhumps are due to modulated 
mass transfer which leads to periodically enhanced dissipation of the stream kinetic energy. 
Smak (2009d) presents evidence for periodically variable irradiation of secondary components 
and proposes that superhumps are due to enhanced dissipation of the stream kinetic energy. 
He outlines specific future modeling efforts that might clarify the overall picture. 



3. Kepler light curve of V344 Lyr 

The bandpass utilized by Kepler is considerably broader than the standard Johnson 
V band. Figure 1 gives a comparison of the two filters |^ Throughout this paper we 
will refer to Johnson V apparent and absolute magnitudes using V and My, and Kepler 
apparent and absolute magnitudes using Kp and Mxp. Figure 2 shows the Kepler light 
curve of V344 Lyr. The observation window contains 19 outbursts showing the sequencing 
SSSLSSSSSSSLSSSSSSS, where S=short outburst and L=long outburst, i.e., a superout- 
burst. Based on photographic observations made by Hoffmeister (1966), Kato (1993) esti- 
mated that the cycle length for normal outbursts was 16±3 d. An estimate of this cycle 
length prior to Kepler was difficult because the normal outbursts of V344 Lyr are short-lived 
and reach only to near Kp = 16 (or V ~ 17). In Table 1 we list the durations (using Kp = 18 
as the cut line) and amplitudes (measured from the local mean quiescence level) of the 19 
observed outbursts (17 normal and two superoutbursts). Both the recurrence times and 
amplitudes of the normal outbursts appear to reach a maximum roughly midway between 
superoutbursts, rather than increasing monotonically between consecutive superoutbursts as 
in the TTIM. We also note that the level of quiescence is ~1 mag higher after the superout- 
burst, and declines with an e— folding time of ~ 50 d. This trend may be associated with 
the cooling of the WD after being heated by the mass accreted during superoutburst (e.g., 
Sion 1995). 



^The Kepler response is taken from 



http: / /keplergo. arc. nasa.gov/ CalibrationResponse.shtml 
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The superoutbursts last ~ 17 d and reach a maximum amplitude /S.Kp — 4.4 — 4.5 
mag. The decline of the superoutbursts are close to linear (plotted as mag versus time, or 
exponential as flux versus time) for ~ 11 d, after which they fall rapidly (~2 d) back to 
quiescence. During this part the decline rate is 0.083 mag d~^, in good agreement with the 
value of 0.094 mag d^^ found by Kato (1993). The amount of time spent at Kp < 16 is ~14 
d, and the duration of the entire plateau is ~12 d. The decay rate is ~12 d mag~^, therefore 
the plateau portion of the decay encompasses ~1 mag — a dynamic range of ~2.5, or slightly 
less than one e— folding (~2.7x). There are also ~0.5 mag shoulders on the superoutburst 
rises lasting ~1 d, which appear to be normal outbursts embedded at the beginning of the 
superoutbursts. 

It is noteworthy that the fast decays of the normal outbursts and superoutburst have 
a faster-than-exponential shape (i.e., concave downward when plotted as magnitude versus 
time) . S84 utilized the Bailey relation between the fast rate of decay in dwarf novae and their 
orbital period to constrain cchot — 0.1 — 0.2. The large errors associated with the AAVSO 
data were consistent with exponential decay. Cannizzo (1994) found that to account for the 
supposedly exponential decays, which correspond to the time during which a cooling front 
traverses the disk, a must vary weakly with radius ( oc r°'"^). The faster-than-exponential 
decays captured by Kepler in V344 Lyr negate this result and indicate that a constant a 
scahng with radius is consistent with the data. 



4. Viscous Plateaus in Dwarf Novae 

The Bailey relation between the fast rate of decay of dwarf nova outbursts and orbital 
period corresponds to the time interval during which the cooling front traverses the disk. 
Some long dwarf nova outbursts have enough dynamic range in V that one can also use the 
slow decay portion, the viscous plateau, to infer an a value. These inferred values agree with 
those determined from the Bailey relation (Warner 1995ab; Cannizzo 2001b, Cannizzo et al. 
2002). 

The decay rates associated with several dwarf novae exhibiting viscous plateaus have 
been fairly well established, and enable one to make a reasonable estimate of the decay rate 
expected in a system with the orbital period of V344 Lyr. The Kepler observation gives a 
decay rate within the superoutburst ~12 d mag~^. We associate this with the viscous time, 
which is basically an e— folding time for a perturbation to the S(r) profile to be damped out. 
For a disk with a S(r) profile approximating steady state, i.e., S(r) oc r~^^^, and for which 
the mass fiow within the disk and onto the WD greatly exceeds the rate of mass addition 
feeding into the outer disk from the secondary star, the viscous time represents roughly a 
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time for the mass of the disk to decrease by a factor e. 

Thus for a given viscous plateau-type outburst, the ratio tpiatcau/^visc represents the 
number of e— foldings by which the mass of the disk decreases due to accretion onto the 
central object. It is instructive to look at examples of viscous decays in several well-studied 
dwarf novae: 

SS Cy^ Long outbursts in SS Cyg (Porb = 6.6 hr) have a duration of ~10 d (Cannizzo & 
Mattel 1992). and the viscous time in the outer disk t^isc ~ 40 d. Thus tpiateau <^ ^visc and 
only about 20% of the disk mass is accreted onto the WD (C93b). 

U Gem: In the 1985 long outburst of U Gem (Porb = 4.25 hr), the decay rate during the 
viscous plateau was ~26 d mag~^ over the ~35 d of the burst (Cannizzo et al. 2002). Since 
^visc ~ ^plateau, ~70% of the disk mass accreted onto the WD. 

V344 Lyr: In the superoutbursts of V344 Lyr (Porb = 2.1 hr) shown in Fig. 2 the locally 
defined decay rate during the plateau remains constant at ~12 d mag~^ over most of the 
~14 d of the burst, thus as with U Gem ~70% of the disk mass accreted onto the WD. 

WZ Sge: In the 2001 superoutburst of WZ Sge (Porb = 81 min), the locally defined decay 
rate during the plateau increased from ~4 to ~12 d mag~^ over the ~20 d of the burst 
(Cannizzo 2001b). The fact that tpiateau ~ Stvisc means that only a few percent of the 
initial disk mass remained at the end of the superoutburst. (An alternative considered 
in the Discussion is that irradiation-induced enhanced mass overflow from the secondary 
star prolonged the superoutburst.) Without a very low acoid — 10^^ there is not enough 
accumulated mass in the disk to account for the outburst's fiuence (Smak 1993) and the 
long recurrence time. 

In summary, one can readily find examples of known systems with viscous type outbursts 
in which the duration of the outburst itself is either greater than, less than, or about the 
same as the viscous time scale evaluated at the outer edge of the accretion disk. The actual 
course taken by a given system depends on the degree of overfilling of the disk in quiescence 
with respect to Mdisk.min- 



^ If one accepts the HST/FGS trigonometric parallax for SS Cyg 6.02 ± 0.46 mas (Harrison et al. 1999) 
which implies D — 166 ± 13 pc, the accretion disk limit cycle model is not able to reproduce dwarf nova 
outbursts for SS Cyg (Schreiber & Lasota 2007, Smak 2010). Harrison et al. (2004) acquire spectra of MK 
standard stars to refine the mean absolute parallax of the reference frame and present a reanalysis of the 
SS Cyg parallax (see their sect. 2.3). They obtain a revised value 6.06 ± 0.44 mas — close to their original 
value. 
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5. Numerical Modeling 

Our numerical model is discussed extensively in C93b and subsequent papers. The basic 
strategy is the time dependent solution of two coupled differential equations, one giving the 
S(r) evolution, and one the evolution of disk midplane temperature T(r). Thus the disk 
is entirely unconstrained as regards deviations from steady state and thermal equilibrium. 
Many previous workers have presented time dependent calculations by solving these equa- 
tions, or similar versions (e.g., S84; Lin, Papaloizou, & Faulkner 1985=L85; Mineshige & 
Osaki 1983=M83, 1985=M85, Mineshige 1986=M86, 1987=M87; Pringle, Verbunt & Wade 
1986=P86; Cannizzo et al. 1986=C86, C93b, 1994, 1998a, 2001a; Angelini & Verbunt 1989; 
Ichikawa & Osaki 1992; Meyer & Meyer-Hofmeister 1984, Ludwig & Meyer 1998; H98, 
Menou, Hameury, & Stehle 1999, Buat-Menard, Hameury, & Lasota 2001, Buat-Menard, & 
Hameury 2002, Schreiber, Hameury, & Lasota 2003, 2004=804). 

Improvements in our code made since the original version include the provision for a 
variable outer disk radius, self-irradiation, and evaporation with various possible laws (e.g., 
Meyer & Meyer-Hofmeister 1994). The diffusion equation governing the evolution of surface 
density is given by (Lightman 1974, Pringle 1981) 



dt r dr 



d_ 
dr 



(^,Sri/2) 



(7) 



where S(r, i) is the surface density and u is the kinematic viscosity coefficient. 



2 a TIT 

We follow the technique of Bath & Pringle (1981) in discretizing the E(r, t) evolution equa- 
tion. 

The thermal energy equation governing the evolution of the midplane disk temperature 
T(r, t) is given by 

dT 2(A-B + C + D) 7^T 1 (9 , , dT 

l^v'^r) -VrjT-, (9) 

at CpE iJ,Cp r or or 

where the viscous heating 



the radiative cooling 



A = -i/Q^s, (10) 
8 



B = aT^, (11) 
3 1 a / „ (9T 
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and 

^ h d f 4acT^dT\ 

D^-—\r — . (13) 

r or \ 3krp or ) 

Of the four terms appearing in the numerator of the first term of thermal energy equa- 
tion, the first two A and B are fairly standard. The term C represents the radial heating 
flux due to turbulent transport, and D is that due to radiative transport. Previous workers 
have used varying forms of C and D. Some early studies adopted C = = (M83), C = 
(L85), or a form of C not expressible as the divergence of a flux (M83, M86, C93b; see H98 
for a critical discussion) . 

The local flow velocity f,., which can vary enormously from inflow to outflow over the 
radial extent of the disk, depending on whether or not transition fronts are present, is given 
by 

where upwind differencing is implemented over the radial grid as follows: 



Sa = Z/i-lSi_irj_i^/^ 

Sb = v-^ir-^'"^ 

Sc = i/i+iSj+irj+i^/^ (15) 



Sb - Sa 



3 Sc — Sb I N 



3 s 



if ( — > and g_ > ) (vj.)j = max(g_,g_|_) 



if ( — > and g_ < 1 {yr)i = min(g_,g+) (17) 



if ( — < and M > 1 1 Mi = q- 

Q+ \<l+\ 

^<0 and ^ 



if ( — < and < 1 ) {vr)i = q+. (18) 
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The scalings for the local maximum and minimum in the E — T relation (from Meyer 
& Meyer-Hofmeister 1981, 1982) are 

^ra^ = 444 g cm-2 r}o'm^°•''a,-°2^ (19) 
where ac-2 — ctcoid/O-Ol, and 

S^i, = 44.4 g cm-2 rlomi'^-^'a^^jl, (20) 
where ah-i — cthot/O.l. The disk midplane temperatures associated with these extrema are 

r(E^a.) = 5275 K (21) 

and 

r(E^,„) = 35900 K rr'mr° °^a;:°jl^ (22) 

The equihbrium temperature scahngs for the stable branches are given in C93b (taken from 
Cannizzo & Wheeler 1984). During transitions when the local disk temperature hes between 
the values associated with the local maximum and minimum in surface density, T'[Einax('")] < 
T{r) < T'[Einin('^)], we use a logarithmic interpolation to obtain the local a value, 

logioQ;(r) = logiottcoid + /, (23) 

where 

^ ^ T(T ■f^Ti'v \ (loglO«hot - logioQicold) ■ (24) 

In tests we also consider a logarithmic interpolation factor / introduced by H98, 

-1 



/ 



1 + 



T 



(logioQihot - logioQicoid), (25) 



as well as a linear interpolation 



a{r) = ttcoid + 7fT7^ — X — ^vv ^ ^"'i^ot ~ "^°i<i)- (26) 

In calculating the effects of the tidal torque from the secondary we follow Smak (1984) 
and later workers who used a formalism developed by Papaloizou & Pringle (1977) in which 
the tidal torque varies as the fifth power of radius. This leads to a depletion in disk material 

= -Cia;— - , (27) 



dt "^njs Va 
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where uj is the orbital angular frequency 27r/Porb) the orbital separation is a, and the specific 
angular momentum jg = (GMir)^/^. Wc do not reset the value of Ci during the course of a 
run as in the TTIM. Ichikawa & Osaki (1994) presented detailed calculations of the strength 
of the tidal torque versus distance at large radii within the disk and found that the power law 
(oc r^) determined by Papaloizou & Pringle is only valid at one radial location in the outer 
disk; further out the torque relation steepens considerably, approaching infinite slope at the 
last non-intersecting orbit. Therefore in some sense the more naive "brick wall" condition 
for the outer edge used by C93b and other workers might be more descriptive. Regardless 
of which prescription is used, if the parameters associated with the outer edge are not 
manipulated during the run, the detailed treatment is not important to the overall model. 
Using a different end point treatment would result in slightly different optimal parameters. 

Given a disk temperature T and surface density E, the disk semithickness h and density 
p are determined from vertical hydrostatic equilibrium to yield 

/.= ^±^2±f^, (28) 

where 

C2 = (29) 
TZT 

C3 = S , (30) 

A* 

the radiation pressure 

Pr = laT\ (31) 

and the density 

P=l (32) 

To determine the computational time step, we take 



At — ft min 



ainSA"^ . fdlnT, 



mm,- — - — , mm, , 

dt ' \ dt 



(33) 



where min^ refers to the minimum over all grid points i, and ft is a small number between 
1/80 and 1/640. The quantities dlnHi/dt and dlnTi/dt are evaluated using the evolution 
equations. 



6. Results 



To calculate the long-term light curves from our numerical trials, we assume local 
Planckian fiux distributions for each annulus, and in each time step integrate over the radial 
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profiles of effective temperature. We compute tlie absolute ii'e]:>/er magnitude Mxp, using tlie 
filter shown in Fig. 1, assuming a face-on disk. To account for the background contribution, 
assumed constant, we add the secondary star and hot spot. For the secondary, we assume 
an M5V star with radius O.221i?0 = 1.54 x 10^° cm and Teg = 2951 K. This contributes 
Mxp = 11.31. For the hot spot, assuming a 12000 K emitting circular region with radius 
10^ cm gives Mxp = 10.41. As a test of our magnitude subroutine, using solar parameters 
Rq = 6.96 X 10^° cm and T^s = 5778 K we find My = 4.85 (the accepted value is 4.83) and 
= 3.83. (The My values for our putative secondary and hot spot are 12.97 and 11.30, 
respectively.) A reasonable inclination i ~ 45° would dim the disk by ~0.4 mag. 

After some experimentation, we found a set of parameters that reproduces approxi- 
mately the overall pattern of supercycle and short outbursts. We take a central O.6M0 
WD accretor, an inner disk radius Tin = 2 x 10^ cm, and a maximum outer disk radius 
'"out = 10^° cm, a mass feeding rate into the outer disk Mt = 2 x 10^^^ Mq yr^^, and a 
values acoid = 0.0025 and Ohot = 0.1. The recurrence time for outbursts varies inversely with 
Mt and Ocoid, and the duration of outbursts varies inversely with ahot (Cannizzo, Shafter 
& Wheeler 1988=CSW, C93b). Also, the ratio of the number of short to long outbursts in 
a long sequence varies in a complicated way on these parameters, and on Tin and rout- The 
strengths of the dependencies of the various time scales — recurrence, burst duration, rise, 
and decay — on the input parameters vary according to location within the multidimensional 
{Mt, acoid, tthot, '"in, '"out, ^i) parameter space (CSW, C93b). 



We use the a interpolation given by eqn. (24). As regards the inner disk radius, the 
value for rjnner is larger than the WD itself; if we take rjnner — ''"wd, we tend to get far 
fewer short outbursts between two successive superoutbursts. This hint of a large inner 
radius may have some observational support in the SU UMa systems (e.g., Howell et al. 
1999). The ratio acoid/«hot = 40 is about 10 times greater than commonly used in dwarf 
nova calculations. The mass transfer rate from the secondary Mt is > 2 times the nominal 
value ~ 0.8 x 10~^^Mq yr~^ found from evolutionary calculations of cataclysmic variables 
below the period gap where the evolution is driven predominantly by gravitational radiation 
(e.g., Kolb, King, & Ritter 1998; Howell, Nelson, & Rappaport 2001). When one takes 
into account the ~20% depletion of the disk mass accompanying each of the short (normal) 
outbursts, however, the effective rate of accumulation of mass in the disk between successive 
superoutbursts decreases from 2 x 10^^*^ M0 yr^^ to ~O.3xlO^^°M0 yr~^. In our models, 
~40— 60% of the disk mass is accreted during a superoutburst. 



In using a complicated numerical model, it is important to carry out testing to gain an 
understanding of systematic effects: 
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6.1. Interpolation between ctcoid and cthot 

Figure 3 shows the effect of using a log versus a Unear interpolation for the a value 
between ctcoid and fthot- The linear interpolation leads to smaller amplitude outbursts with 
little or no quiescent intervals. The decays are also much slower, and concave upwards when 
plotted as mag versus time (i.e., with a functional form that is slower-than-exponential) . 
Figure 4 shows the effect of using the H98 interpolation. H98 adopted the constant To = 
25000 K in their scaling; we also look at smaller values 18750 K and 12500 K. The smaller 
To values produce superoutbursts which have the recognizable plateau and subsequent fast 
decay, but for the Tq = 12500 K run the fast decays are slower-than-exponential. Hameury 
(2002) compared the two cases Tq = 24000 K and To = 8000 K and found httle difference 
between the resultant light curves (see his Fig. 3) in his model. 

Figure 5 depicts the different interpolations as well as parameterizations of the fractional 
ionization ^. The MRI is thought to mediate the strength of a through turbulent transport, 
and the effect of shearing on the weak magnetic field embedded in the gas is more efficient 
as ^ increases. The H98 scaling for To = 12500 K does a reasonable job in approximating 
the p = 10^^ g cm~^ curve for ^, and may be more physical than the other scalings. 

6.2. Number of Grid Points N 

Early time dependent accretion disk limit cycle studies were concerned mainly with 
showing the ability of the disk, under the conditions of an imposed limit cycle in each annulus, 
to undergo collective, global oscillations in which a major part of the disk participated. 
The total numbers of grid points N was small and the comparison with observational data 
minimal (e.g., N ^ 20 - 25: S84; N = 35: L85; N ^ 22 - 44: M85; N = 17: M86; 
N — 40: P86; N — 44: C86). One can see obvious irregularities in many of the early light 
curves which are characteristic of numerical instability. C93b undertook the first detailed 
investigation regarding the number of grid points, and found that at least several hundred 
were required to obtain reliable light curves. H98 utilized both uniform and adaptive mesh 
codes to study the properties of the transition fronts, and found a requirement oi N > 800 
for a uniform-equivalent mesh to resolve the fronts adequately. 

Figure 6 shows a detail of one superoutburst for models with — 200 to 1600 grid 
points. In term of the long term light curves, even with > 1000 the sequencing of long 
and short outbursts is not perfectly stable using our explicit code. An implicit code such 
as that described in H98 may be more numerically stable. In terms of individual outburst 
profiles, the rise times we calculate for both the normal and superoutbursts are too fast. In 
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addition, our outbursts are too bright by ~0.5 mag. A nominal inclination i ~ 45° would 
make our disks dimmer by ~0.4 mag, bringing them close to the level observed in outburst. 
For comparison we also plot the V344 Lyr data, where D = 620 pc was assumed in converting 
to absolute magnitude (found by Ak et al. 2008 using a period-luminosity relation). There 
is a slight increase in superoutburst duration from = 200 to A^ = 800 due to the fact 
that more mass accumulates as the supercycle lengthens. In addition, the shoulder in the 
light curve also becomes longer. For the highest A^ run the duration of the superoutburst is 
still shorter than that seen in V344 Lyr, while the shoulder is too long. The superoutburst 
calculated using the H98 scaling with Tq = 12500 K better reproduces the observed duration 
of the V344 Lyr superoutburst, and also has a shorter shoulder with less structure than the 
other calculations. However, it suffers from a rise time that is too short, and the fast decay 
is now slower-than-exponential. Also, the viscous decay rate, which has been optimized to 
the observed value using the eqn. (24) scaling, is now too fast. 



6.3. Thermal Energy Equation 



Considering the many forms of the thermal energy equation used by previous workers, 
and the potential limitations of employing a one-dimensional hydrodynamical formalism, we 
now look at the effects of switching off various terms in the thermal energy equation in an 
attempt to ascertain the relative importance of each term. Figure 7 shows the effect of setting 
equal to zero: (i) C, (ii) D, (iii) C and D, (iv) the first advective term in eqn. and (v) 
the second advective term in eqn. (g. In the C = run one sees a rounded superoutburst 
and only three short outbursts between two superoutbursts. The (radiative) D term is less 
influential, as the superoutburst profile is relatively unchanged from the first panel. Setting 
C = D = produces a single dip in the rise to superoutburst, and a longer viscous plateau 
(i.e., exponential decay). The first advective term appears to be of comparably minimal 
importance as the D term; i.e., the omission of either has a negligible effect. The second 
advective term is the most influential of all terms examined: without it one has much longer 
superoutbursts — about 70% of the accumulated disk mass is accreted, which is about the 
same as in V344 Lyr. However, the viscous plateau still comprises less than half of the 
superoutburst, rather than the > 90% seen in V344 Lyr. 



6.4. Detuning of Optimal Parameters 



There is a complicated interplay between the different parameters that enter into the 
overall results. Figure 8 shows the effect of varying the inner disk radius. For smaller inner 
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disk radii, one has a larger djTiamic range in router /''^inner and as a result there tends to be a 
more complete emptying of the disk during outbursts, and therefore less of a tendency to have 
a series of short outbursts leading up to a superoutburst. The net interaction among all the 
model parameters is such that the observed sequencing for V344 Lyr is roughly reproduced 
with Tinner — 2 X 10^ cm. Figure 9 shows the effect of varying the outer disk radius. For 
larger router values not only does the number of short outbursts increase, but also the shape 
of the superoutburst distorts significantly from the viscous plateau shape. 



6.5. Relevance of Results to the V344 Lyr light curve 

No one light curve we have shown reproduces all features of the observed light curve 
shown in Fig. 2. For instance, there is an atypically long quiescent period at the beginning 
of the observed light curve. The recurrence time for outbursts varies inversely to varying 
degrees with ttcoid, «hot, and Mt (CSW). Therefore if any one of these parameters were 
anomalously small, that could delay an outburst. Also, the pre-superoutburst normal out- 
bursts are brighter than the post-superoutburst normal outbursts, and the pre-superoutburst 
quiescent level is fainter than the post-superoutburst quiescent level. As noted earlier, the 
disk instabihty model is ill-equipped to make any statement about the quiescent state. How- 
ever, one possibility might be that were elevated for a while after the superoutburst; 
that could lead to crowding of the short outbursts, and an increase in the quiescent level. 
There may be other causal relations among acoid, «hot, and Mt, which are important. 

An analysis of the phase variations in the quiescent-state oscillations in the light curve 
at the beginning of Fig. 2 light curve show these are likely to be negative superhumps. As 
such they indicate accretion onto an accretion disk that is tilted a few degrees out of the 
orbital plane and precessing slowly in the retrograde direction (Barrett et al. 1988; Foulkes 
et al. 2006; Wood & Burke 2007, Wood et al. 2000, 2009). Instead of impacting the rim 
of the disk, the accretion spot sweeps across first one face of the tilted disk and then the 
other during one orbit, impacting the rim only twice during one cycle, along the line of 
nodes. Thus, the accretion mass and energy is deposited throughout a range of radii in 
the disk, and this may act to stabilize the disk against DN outbursts. The slow rise of the 
first dwarf nova oscillations after the long quiescent state indicates that this is an inside-out 
outburst, consistent with the deposition of mass primarily in the inner disk. Clearly these 
data motivate further study of the thermal-viscous limit cycle instability in tilted disks. 



The primary interpolation scheme we considered, given in eqn. (24), has clear problems 
in making good superoutbursts. The model superoutbursts have shoulders which comprise 
about a third to half of the total viscous plateau, rather than the ~10% observed. However, 
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for the H98 scaling with Tq = 12500 K, the shoulders are appropriately short, ~2 d, and the 
viscous plateaus much longer. The rate of decay of the viscous plateau is too fast, but that 
is because the input parameters have been optimized for the scaling given in eqn. (24). 



7. Discussion 

In contrast to the thermal-tidal instability of Osaki (Ichikawa & Osaki 1992, 1994; 
Ichikawa, Hirose, & Osaki 1993) in which the tidal torque is artificially enhanced by a factor 
of 20 when the outer disk expands beyond the point of 3:1 commensurability with the orbital 
period, our detailed modeling shows that the regular and superoutbursts can produced with 
a constant set of parameters. In retrospect, one may legitimately question the motivation 
for the increased tidal torque accompanying the expansion of the disk, considering that the 
ultimate driving force behind the tidal torque is the gravitational field of the secondary, 
which remains unchanged. We share a common viewpoint with van Paradijs (1983) in that 
we regard the superoutbursts in the SU UMa systems as an extension of the long outbursts 
in DNe above the period gap, such as U Gem and SS Cyg (Cannizzo & Mattel 1992, 1998; 
C93b). Continued monitoring of V344 Lyr should reveal the stability of the long term 
sequencing, i.e., variations in the number of shorts sandwiched between two longs. 

The shoulders on the rise to superoutburst in the models have more structure than the 
shoulder in the observed superoutburst. The shoulders arise because the initial tendency 
of the triggering of thermal instability is to produce a short outburst. However, due to the 
long-term net build-up of mass in the disk at large radii between successive short outbursts, 
there is actually enough material accumulated to support a superoutburst. It takes a little 
longer however for the outward moving heating front to progress to larger radii and to fully 
access this additional store of matter, which manifests itself as a ring of enhanced surface 
density in the outer disk. This explains the fact that many observed superoutbursts appear 
to start as normal (short) outbursts. The subsequent rise after the shoulder is due to the 
full radial extent of the disk out to the outer edge making a transition to the hot state. In 
contrast, for the short (normal) outbursts, not enough disk mass has accumulated at larger 
radii by the time of the outburst to support the heating transition front propagating to large 
radii, therefore the heating front stalls at some r < router and is reflected as a cooling front, 
leading to pointy maxima rather than viscous plateaus. 

Our testing has delineated a range of allowable parameters: 

(1) The logarithmic interpolation for a produces well-separated outbursts, with concave 
downward decays, as observed. The linear interpolation produces outbursts with minimal 
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quiescence and concave upward decays. Given the physical underpinnings of the MRI, the 
logarithmic interpolation is better motivated since the MRI strength should depend on the 
partial ionization fraction, which is a strongly increasing function of temperature. 

In terms of the superout burst, the H98 scahng with Tq — 12500 K produces the best 
overall results. The shoulder feature only persists ~ 2 d, as observed, and the viscous plateau 
is longer. About 70% of the stored mass is accreted. The main deficiencies are that the fast 
decay (i.e., that following the viscous plateau) is slower-than-exponential, and the rise times 
are too fast (as with the other calculations). 

(2) The minimum number of grid points necessary is > 10^ (in accord with C93b and 
H98). Even for high A^, the sequencing is not perfectly stable. 

(3) Of the terms entering into the thermal energy equation, the radial viscous transport 
term C and the second advective term, Vr{dT/dr) appear to be the most important. The 
radial radiative fiux term D and the first advective term, involving d{rVr)/dr, are of lesser 
importance. 

(4) The values for acoid = 0.0025, cthot = 0.1, and Mr = 2 x yr~^ are mandated 
by the observed spacing and duration of short outbursts and the superout burst. Increasing 
their values decreases the outburst time scales. The inferred values for cicoid and ahot are 
to be viewed merely as best-fit parameters for V344 Lyr; they do not allow us to draw any 
deep conclusions as regards other systems. 

(5) The values Tinner and router strongly mediate the number of short outbursts between 
two superout bursts. For N = 800 and ft = 1/160 we have ~5 shorts between two longs for 
'"inner = 2 X 10^ cm, or ~3i?wD fo^ MyvD = O.GMq. Observations by Howell et al. (1999) 
seem to favor an evacuated inner disk. Several possibilities are discussed by these workers, 
including the possibifity of weak intermediate polars in the SU UMa systems. The outbursts 
in our models are triggered at intermediate disk radii, therefore the outer disk radius has 
minimal effect on the outburst recurrence times and durations. Increasing router from 1 x 10^° 
cm to 2 X 10^° cm increases the number of shorts between to successive longs, and also distorts 
the profile of the superout burst. 

With tthot = 0.1, we match the rate of decay of the viscous portion of the superout burst, 
and also find a ~0.5 mag shoulder on the rise to superout burst. It is not clear why previous 
detailed time dependent modelers did not find this shoulder (e.g., S04). The main deficiency 
in the models with the standard logarithmic acoid cthot interpolation is their failure to 
produce a sufficiently long viscous plateau on the superout burst. The observed plateau has 
a decay rate of ~12 d mag^^ and lasts for ~12 d, which corresponds to ~1 e— folding in the 
decay of the disk mass, meaning about 60% of the disk mass should be accreted. Our viscous 
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plateau lasts ~6— 8 d, with proportionally less mass accreted. However, the light curve using 
the H98 scaling with Tq = 12500 K (which appears to be physically motivated insofar as it 
tracks ^) has substantially longer superout bursts, ~ 10 d, with appropriately short shoulders 
as observed. The reason for the shorter shoulders is that the temperature within the outward 
moving heating front accompanying superoutburst onset is between about 10^ and 2 x 10"' 
K. As can be seen in Figure 5 by comparing the left-most green curve (i.e., the H98_12500 
scahng) with the "log" curve, for the temperature range 10^ K < T < 2 x 10^ K, a is larger 
in the H98_12500 scaling and therefore the heating front travels faster since its speed varies 
with a (Menou, Hameury, & Stehle 1999). The duration of the shoulder feature corresponds 
to the interval during which the heating front is moving outward and finding that enough 
mass is present at large radii to support a superoutburst. This causes the dip and subsequent 
rise to the plateau of the superoutburst. 

The values for ctcoid, cchot, and Mt have already been optimized to get the supercycle 
and normal cycle periods of ~100 d and ~10 d, respectively. One suspects that the scalings 
for Sinax(«coid) and Smin(tthot) Sire too simplistic and fail to take into account the actual disk 
physics, particularly that of the quiescent state, which is presently unknown. One could 
envision introducing an ad hoc multiplicative parameter for Smax, for instance to artificially 
enhance it, but then one would have to re- adjust all the other parameters, in particular 
increasing Mt which is already uncomfortably large. In view of these considerations we 
consider our experiment moderately successful in that we are able to reproduce the duration 
of the viscous plateau within a factor of two in an over-constrained system for which ctcoid 
and cchot have already been determined by the outburst properties. 

What about having enhanced mass transfer from the secondary prolong the duration 
of the superoutburst? Smak (2004) presented observational evidence for enhanced hot spot 
brightness near superoutburst maxima which he interpreted as being due to enhanced mass 
transfer brought about by the strong irradiation of the secondary. Osaki & Meyer (2003) pre- 
sented a theoretical study of the effects of irradiation on the photosphere of the secondary star 
in SU UMa systems and concluded that mass transfer could not be enhanced by a significant 
amount because the Coriolis force prevents the formation of a circulation fiow transporting 
heat toward LI. Viallet & Hameury (2007= V07, 2008= V08) present two-dimensional, time 
dependent calculations of the surface flow of material of irradiated secondaries and flnd the 
physical situation to be signiflcantly different from that envisioned by either Osaki & Meyer 
or Smak. The Coriohs force leads to the formation of a circulation flow from the high lat- 
itude regions to the vicinity of the LI point, but rapid coohng of the gas as it enters the 
shielded equatorial region prevents significant heat transport (V07). V07 conclude that any 
resulting mass transfer rate enhancement is likely to be moderate and unable to account for 
the duration of long outbursts in dwarf novae. V08 look at two effects that could potentially 
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induce secondary mass transfer rate variations in dwarf novae and soft X-ray transients: ir- 
radiation of the secondary by the hot outer disk rim during outburst, and scattered radiation 
by optically thin outflowing material. They conclude that for dwarf nova parameters the 
effects on enhanced mass transfer are negligible. For completeness we note that some of the 
assumptions of V07 and V08 may not apply if the disk is warped or tilted, in which case the 
LI point could in principle be directly irradiated. 

S04 presented detailed time dependent calculations of the accretion disk evolution in 
both the TTIM and the EMTM, and in the end concluded that the EMTM was favored. 
Their arguments arc rather general, however, and they state "... we have not proven the 
EMTM to be correct nor the TTIM not to work." Both models have a number of free 
parameters in addition to those already present in the standard limit cycle model. For the 
TTIM one has the time delay between the attainment of the 3:1 radius by the outer disk edge 
and the onset of enhanced tidal torque, in addition to the functional form and amplitude 
associated with the increase. Similarly the EMTM has as free parameters the degree of 
augmentation of Mt during superout burst, as well as its functional form with time. 

In addition, there may be a problem with the viscous decay in the EMTM: The V344 
Lyr superout burst maintained a closely exponential decay over ~1 mag, whereas our models 
as well as others reveal a definite slower-than-exponential trend (i.e., concave upward when 
plotted as magnitude versus time) in viscous plateaus covering that much dynamic range. 
Although the viscous decays can be lengthened in the EMTM, they tend to have noticeable 
deviations from exponcntiality, in contrast to the V344 superoutburst. The fact that the 
short orbital period system WZ Sge did show a deviation from exponcntiality in its 2001 
outburst, which also spanned a greater dynamic range than in the V344 Lyr superoutburst 
(~ 2 mag versus ~1 mag), may indicate that irradiation-induced enhanced mass transfer is 
a factor at much shorter orbital periods than for V344 Lyr. In addition, the amphtude of 
the ~0.5 mag shoulder on the rise to superoutburst is roughly matched in our calculations; 
with enhanced mass transfer the shoulder amplitude is too great (> 1 mag). We consider it 
fortunate that we have matched the decay as well as we have since the scalings for the extrema 
in the limit cycle model, in particular Sinax(ctcoid)) ai'e most probably gross simplifications 
to a more complicated situation that is currently beyond our understanding. 



8. Conclusion 

We have presented detailed time dependent calculations of the accretion disk limit cycle 
mechanism with application to the Kepler light curve of V344 Lyr. We find that the basic 
disk instability model is able to account for the mixture of normal and superoutbursts with 
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a minimal amount of fine-tuning. Given the additional free parameters associated with 
both the thermal-tidal model and the enhanced mass transfer models, we do not attempt 
to examine their possible ramifications in this work. It may be that the plain vanilla disk 
instability acting alone is sufficient to explain most aspects of the long term V344 Lyr light 
curve. We also note that the normal outbursts in V344 Lyr do not show a monotonic 
increase in recurrence time and amplitude between successive superoutbursts as predicted 
by the TTIM (see Ichikawa, Hirose, & Osaki 1993, Fig. 1; Osaki 2005, Fig. 3), but rather a 
maximum in these quantities near the middle of a supercycle. The models we present also 
do not show the observed variations. 

One can in some sense view the alternation of many successive short outbursts with 
a sup erout burst as an extensive of the pattern of short and long outbursts seen in the 
longer orbital period systems such as U Gem and SS Cyg (van Paradijs 1983); for SU 
UMa parameters one has many more short outbursts between two longs than in a system 
like SS Cyg, which exhibits LSLSLS... most of the time (Cannizzo & Mattel 1992). For 
our best-fit model parameters, we find that > 800 grid points are required for numerical 
stability, with a resulting pattern in which 5—7 short outbursts occur between two successive 
superoutbursts (for ft = 1/640). Overall, our outbursts are at about the same brightness 
as those observed in V344 Lyr for an intermediately inclined disk. Having constrained 
CKcoid s-nd tthot so as to produce the overall recurrence time scales of short outbursts and 
superoutbursts, there are several obvious deficiencies: the rise times of all outbursts are 
too fast, the superoutburst duration is too short, and the duration of the shoulder is too 
long and has an extra scalloped-out feature prior to the observed shoulder. However, the 
more physically motivated H98 scaling with Tq = 12500 K does a better job on making 
superoutbursts with minimal shoulders as observed, lasting ~2 d, and also longer viscous 
plateaus, ~ 10 d vs. the observed ~ 12 d. The normal outburst embedded at the start of a 
superoutburst is due to the initial triggering of the thermal instability and outward movement 
of the heating front; the subsequent dip and rise to superoutburst maximum comes about 
because the heating front continues to move outward, and subsequently to access the stored 
material at large radii in the disk which accumulates where the heating front spike stalls at 
the end of a normal outburst. 

Even as we struggle to account for the finer points of the light curve, such as the shoulder 
on the superoutburst, we feel compelled to step back briefly from our myopic scrutiny of the 
details and give some modicum of thanks and appreciation for the fact that the Kepler data 
represent such a huge advance over what was available even as recently as one year ago. For 
long term dwarf nova light curves, the only data sets that had been available were those 
amassed by the AAVSO and similar organizations. While these data have certainly been 
useful and yielded many interesting results, one could only study in detail the brightest 
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systems, like SS Cyg and U Gem. In addition, the data quality was uneven and the time 
sampling coarse. When new data become available that arc an order of magnitude better 
than the previous data, it becomes possible to make discoveries that are different in "kind" , 
not just in "degree". The Kepler data will provide a physical touchstone for modelers with 
the potential to substantially improve our understanding of the physics of accretion disks; 
our work represents only the first small step in this regard. 

We thank James Bubeck (under the SESDA 11 contract) at Goddard Space Flight Center 
for assistance in providing local computational resources, and Marcus Hohlmann and Patrick 
Ford from the FIT High Energy Physics group and the Domestic Nuclear Detection Office in 
the Dept. of Homeland Security for making additional resources on a Linux cluster available. 

We also thank the anonymous referee whose comments helped improve and clarify the 
paper. We acknowledge the contributions of the entire Kepler team. 
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Table 1: Duration and Amplitude of Outbursts 
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Figure Captions 



Fig. 1. — A comparison of tlie Jotinson V band filter (red) with the Kepler filter (blue). 

Fig. 2. — Kepler light curve of V344 Lyr over 270 d showing outbursts every ~10 d, with 
two superoutbursts. The flux is measured in e~ cadence"^, where the cadence (time between 
integrations) is 1 m. The full time resolution is shown. To obtain visual magnitude, a 
correspondence of Kp = 12 to 10^ e~ cadence"^ was adopted. (This conversion assumes an 
unreddened G2 star.) 

Fig. 3. — The effect of using a logarithmic interpolation (eqn. [24]) versus a linear one 
(eqn. [26]) for a. For these runs = 800 and ft = 1/160. Shown are the light curves for 
the logarithmic case {top panel) and the linear case {second panel), the accompanying disk 
masses {third and fourth panels), and the radius of the hot/cold interface rHc,io during times 
when a transition front is present {fifth and sixth panels). For comparison the light curves 
for the disk-only contribution (i.e., excluding secondary and hot spot) are shown as dotted 
curves in the first and second panels. The disk masses are in units of 10^^ g and rHc,io is in 
units of 10^° cm. 

Fig. 4. — The effect of using the H98 logarithmic factor / (eqn. [25]). For these runs 
= 800 and ft = 1/160. Shown are the light curves for Tq = 25000 K {top panel), 18750 
K {second panel), and 12500 K {third panel), the accompanying disk masses {fourth through 
sixth panels), and the radius of the hot/cold interface rncio during times when a transition 
front is present {seventh through ninth panels). The disk masses are in units of 10^^ g and 
'^HCio is in units of 10^° cm. 

Fig. 5. — Shown are the different acoid — tthot interpolations. The blue lines show the scalings 
for T(Sjnax) and T(Sjnin)- The black curves indicate the hnear (eqn. [26]) and logarithmic 
(eqn. [24]) interpolations. The green curves show the H98 scaling (eqn. [25]), where the 
three curves (left to right) represent Tq = 25000 K, 18750 K, and 12500 K. The three red 
curves show a measure of the partial ionization; they are the relative contribution of electron 
pressure to total pressure as shown in Fig. 16 of C93b, for p{g cm~^)= 10~^, 10"'', and 
10~^ (left to right). The values given in C93b have been multiplied by 2ahot to take into 
account the fact that (1) Pe/P 0.5 in the limit of complete ionization, and (2) in the high 
temperature limit a — j- ahot = 0.1. 

Fig. 6. — Model light curves showing the detail of a superoutburst for separate calculations 
in which A^ = 200 {short dashed), 400 {long dashed), 800 {short dash-dotted), and 1600 {long 
dash-dotted). The green curve shows a superoutburst from the Tq = 12500 K H98 scaling 
light curve of Fig. 4. For comparison the averaged Kepler data for the superoutbursts in 



Fig. 2 is shown {first: blue, second: magenta), where the original At — 1 min cadence has 
been block-averaged to At = 4 hr and a distance of 620 pc was assumed in converting from 
Kp to Mxp. The two slanted red lines indicate a decay slope of 12 d mag~^. The second 
panel shows the disk mass (in 10^^ g). 

Fig. 7. — The effect of turning off different terms in the energy equation. For these runs 
N — 800 and ft — 1/160. Shown are light curves using: the full equation {top panel), the 
turbulent transport C = {second panel), the radiative transport D = {third panel), both 
C = and D = {fourth panel), the first advective term set to zero {fifth panel), and the 
second advective term set to zero {sixth panel). 

Fig. 8. — The effect of the inner disk radius. For these runs N — 800 and ft — 1/160. The 
four panels show the resultant hght curves for Tinner — 2 x 10^ cm {top panel), 1.5 x 10^ cm 
{second panel), 1 x 10^ cm {third panel), and 0.5 x 10^ cm {fourth panel). 

Fig. 9.— The effect of the outer disk radius. For these runs = 800 and ft = 1/160. The 
six panels show the resultant hght curves for Tinner = 0.8 x 10^° cm {top panel), 1.0 x 10^° 
cm {second panel), 1.25 x 10-*^° cm {third panel), 1.5 x 10^*^ cm {fourth panel), 1.75 x 10^° cm 
{fifth panel), and 2.0 x 10^° cm {siaih panel). 
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